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Abstract 

We present a new measurement of the B + meson differential cross section da/dpx at y/s = 1960 
GeV. The data correspond to an integrated luminosity of 739 pb _1 collected with the upgraded 
CDF detector (CDF II) at the Fermilab Tevatron collider. B + candidates are reconstructed through 
the decay B + — ► J/ip K + , with J/ip — ► u + . The integrated cross section for producing B + 
mesons with px > 6 GeV/c and \y\ < 1 is measured to be 2.78 ± 0.24 fib. 

PACS numbers: 13.85.Ni, 12.38.Qk, 13.25.Hw, 14.40.Nd 
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I. INTRODUCTION 



Measurements of the bottom quark production cross section at the Tevatron collider 
probe the ability of perturbative QCD to predict absolute rates in hadronic collisions. At 
the perturbative level, calculations of the hard scattering cross sections have been carried 
out at next-to-leading order (NLO) [1] and also implemented with logarithmic p^/mi, 1 
corrections evaluated to next-to- leading logarithmic accuracy (NLL) [2]. In both cases, 
these QCD predictions are affected by large theoretical uncertainties such as the dependence 
on the choice of the renormalization and factorization scales and the 6-quark mass [3, 4]. 
Accurate measurements could help in improving the theoretical prediction. Unfortunately, 
as noted in Ref . [5] , measurements of the 6-quark cross section at the Tevatron appear to be 
inconsistent among themselves. Reference [5] uses the prediction of a NLO calculation [1] 
implemented with a non-perturbative model for the 6-quark fragmentation 2 in order to 
compare all measurements performed at the Tevatron. This calculation predicts cr{p B+ > 
6 GeV/c, \y\ B+ < 1) = 0.9 fib. Previous measurements [9, 10] performed by the CDF 
collaboration at y/s = 1.8 TeV yield a{j>% + > 6 GeV/c, \y\ B+ < 1) = 2.66±0.61 and 3.6±0.6 
yub, respectively. The ratios of these measurements to the NLO prediction are (2.9 ± 0.67) 
and (4.0 ± 0.6), respectively. In contrast, the ratios of the CDF and DO measurements of 
the b cross section, that are not based upon the detection of J/ip mesons [11-15], to the same 
theoretical prediction have an appreciably smaller average (2.2 with a 0.2 RMS deviation [5]). 
The cause of the inconsistency could be experimental difficulties inherent to each result or 
some underlying, and not yet appreciated, production of new physics. Therefore, it is of 
interest to clarify the experimental situation. 

This paper presents a new measurement of the B + production cross section that uses fully 
reconstructed — > J/ip decays. We follow closely the experimental procedure used 
in Refs. [9, 10], but we simplify the analysis selection criteria in order to reduce systematic 
uncertainties. The B + production cross section is the ratio of the number of observed B + 
candidates to the product of the detector acceptance, integrated luminosity, and branching 

1 Mass (m,b) and transverse momentum (p^) of the bottom quarks involved in the hard scattering. 

2 This calculation uses a 6-quark mass of rrif, = 4.75 GeV/c 2 , renormalization and factorization scales 
/in = jip = \Jvt + m b> the MRSD [6] fit to the parton distribution functions (PDF), and a fragmentation 
fraction /„ = 0.375. The fragmentation model is based on the Peterson fragmentation function [7] with 
the e parameter set to 0.006 according to fits to e + e~ data [8]. 
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fraction of the decay B + — > J/ip K + with J/ip — > /i + //~. Section II describes the detector 
systems relevant to this analysis. The data collection, event selection, and B ± reconstruction 
are described in Sec. III. In Sec. IV, we evaluate the detector acceptance and derive the 
total and differential B + cross section. Our conclusions are presented in Sec. V. 

II. THE CDF II DETECTOR 

CDF is a multipurpose detector, equipped with a charged particle spectrometer and a 
finely segmented calorimeter. In this section, we describe the detector components that are 
relevant to this analysis. The description of these subsystems can be found in Refs. [16-22]. 
Two devices inside the 1.4 T solenoid are used for measuring the momentum of charged 
particles: the silicon vertex detector (SVX II) and the central tracking chamber (COT). 
The SVX II consists of double-sided microstrip sensors arranged in five cylindrical shells 
with radii between 2.5 and 10.6 cm. The detector is divided into three contiguous five-layer 
sections along the beam direction for a total z coverage 3 of 90 cm. The COT is a cylindrical 
drift chamber containing 96 sense wire layers grouped into eight alternating superlayers of 
axial and stereo wires. Its active volume covers \z\ < 155 cm and 40 to 140 cm in radius. The 
central muon detector (CMU) is located around the central electromagnetic and hadronic 
calorimeters, which have a thickness of 5.5 interaction lengths at normal incidence. 

The CMU detector covers a nominal pseudorapidity range < 0.63 relative to the 
center of the detector, and is segmented into two barrels of 24 modules, each covering 15° 
in 0. Every module is further segmented into three submodules, each covering 4.2° in and 
consisting of four layers of drift chambers. The smallest drift unit, called a stack, covers 
a 1.2° angle in <fi. Adjacent pairs of stacks are combined together into a tower. A track 
segment (hits in two out of four layers of a stack) detected in a tower is referred to as a 
CMU stub. A second set of muon drift chambers (CMP) is located behind an additional 
steel absorber of 3.3 interaction lengths. Muons which produce a stub in both CMU and 
CMP systems are called CMUP muons. 

3 In the CDF coordinate system, 9 and <j> are the polar and azimuthal angles of a track, respectively, defined 
with respect to the proton beam direction, z. The pseudorapidity r\ is defined as — log tan(0/2). The 
transverse momentum of a particle is pt — P sin(0). The rapidity is defined as y — l/2-\og((E+p z )/(E — 
Pz j), where E and p z are the energy and longitudinal momentum of the particle associated with the track. 
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The luminosity is measured using gaseous Cherenkov counters (CLC) that monitor the 
rate of inelastic pp collisions. The inelastic pp cross section at y/s = 1960 GeV is scaled 
from measurements at y/s = 1800 GeV using the calculations in Ref. [23]. The integrated 
luminosity is determined with a 6% systematic accuracy [24]. 

CDF uses a three-level trigger system. At Level 1 (LI), data from every beam crossing 
are stored in a pipeline capable of buffering data from 42 beam crossings. The LI trigger 
either rejects events or copies them into one of the four Level 2 (L2) buffers. Events that pass 
the LI and L2 selection criteria are sent to the Level 3 (L3) trigger, a cluster of computers 
running speed-optimized reconstruction code. 

For this study, we select events with two muon candidates identified by the LI and L2 
triggers. The LI trigger uses tracks with px > 1.5 GeV/c found by a fast track processor 
(XFT). The XFT examines COT hits from four axial superlayers and provides r — <fi informa- 
tion. The XFT finds tracks with p T > 1.5 GeV/c in azimuthal sections of 1.25°. The XFT 
passes the tracks to a set of extrapolation units that determines the CMU towers in which 
a CMU stub should be found if the track is a muon. If a stub is found, a LI CMU primitive 
is generated. The LI dimuon trigger requires at least two CMU primitives, separated by at 
least two CMU towers. At LI, there is no requirement that muons have opposite charge. 
During the data-taking period in which the dimuon sample used for this analysis was col- 
lected, the Tevatron luminosity has increased from 1 to 100 xlO 30 cm" 2 s -1 . Accordingly, the 
L2 trigger, that started with no additional requirement, has incrementally required dimuons 
with opposite charge, opening azimuthal angle 5<p < 120°, and pr > 2 GeV/c. All these 
trigger requirements are mimicked by the detector simulation on a run-by-run basis. At L3, 
muons are required to have opposite charge, invariant mass in the window 2.7 — 4.0 GeV/c 2 , 
and \5z \ < 5 cm, where z is the z coordinate of the muon track at its point of closest 
approach to the beam line in the r — plane. These requirements define the J/ip — > yU + /i" 
trigger. 

We use two additional triggers in order to verify the detector simulation. The first 
trigger (CMUPpH) selects events with at least one LI and one L2 CMUP primitive with 
Pt > 4 GeV/c, and an additional muon found by the L3 algorithms. Events collected with 
this trigger are used to measure the muon trigger efficiency The second trigger (//— SVT) 
requires a LI CMUP primitive with p? > 4 GeV/c accompanied by a L2 requirement of 
an additional XFT track with p T > 2 GeV/c and displaced from the interaction point. 



10 



These events are used to verify the muon detector acceptance and the muon reconstruction 
efficiency. 

III. DATA SELECTION AND B ± RECONSTRUCTION 

We search for B ± — > J / ipK ± candidates in the data set selected by the J/ip — > fi + fi~ trig- 
ger. Events are reconstructed off-line taking advantage of more refined calibration constants 
and reconstruction algorithms. 

The transverse momentum resolution of tracks reconstructed using COT hits is 
°(Pt)/Pt - 0.0017 [GeV/c]- 1 . COT tracks are extrapolated into the SVX II detector 
and refitted adding hits consistent with the track extrapolation. Stubs reconstructed in the 
CMU detector are matched to tracks with pt > 1.3 GeV/c. A track is identified as a CMU 
muon if Ar0, the distance in the r — <fi plane between the track projected to the CMU 
chambers and a CMU stub, is less than 30 cm. We also require that muon-candidate stubs 
correspond to a LI CMU primitive, and correct the muon momentum for energy losses in 
the detector. 

We search for J/ip candidates by using pairs of CMU muons with opposite charge, and 
Pt > 2 GeV/c (this requirement avoids the region of rapidly changing efficiency around 
the trigger threshold). The invariant mass of a muon pair is evaluated by constraining 
the two muon tracks to originate from a common point in three-dimensional space (vertex 
constraint) in order to improve the mass resolution. All muon pairs with invariant mass in 
the range 3.05 — 3.15 GeV/c 2 are considered to be J/ip candidates. 

If a J/ip candidate is found, we search for B ± mesons by considering all remaining charged 
particle tracks in the event as possible kaon candidates. As in previous measurements [9, 10], 
we select tracks with px > 1-25 GeV/ c and with | <5^o I < 1-5 cm with respect to the Zq position 
of the J/ip candidate. We require that kaon-candidate tracks have at least 10 hits in both 
COT axial and stereo superlayers. This limits the pseudorapidity acceptance to \q\ < 1.3. 
The invariant mass of the ji + n~ K ± system is evaluated constraining the corresponding 
tracks to have a common origin while the /z + \i~ invariant mass is constrained to the value of 
3.0969 GeV/c 2 [25]. As in Refs. [9, 10], we select B ± candidates with p T > 6 GeV/c. From 
the pseudorapidity acceptance of CMU muons < 0.8) and the Pt cuts on the /i ± and 
B ± transverse momenta, it follows that: (1) no kaon from B ± decays is emitted at \r\\ > 1.3; 
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(2) the reconstructed B ± candidates have rapidity \y\ < 1. 

In contrast with the analyses in Refs. [9, 10], we do not require the proper decay length 
of the B ± candidates to be larger than 100 pim. By doing so, we avoid two large sources 
of systematic uncertainty: (1) the simulated efficiency of the SVX II detector; (2) the 
dependence of the decay length distribution on the simulated SVX II resolution and B ± 
transverse momentum distribution. The invariant mass distribution of all B ± candidates 
found in this study is shown in Fig. 1. 

IV. DIFFERENTIAL CROSS SECTION 

To measure the differential cross section, we divide the sample of B ± candidates into 
five p T bins: 6 - 9, 9 - 12, 12 - 15, 15 - 25, and > 25 GeV/c [9, 10]. In each p T bin, we 
fit the invariant mass distribution of the B^ candidates with a binned maximum-likelihood 
method to determine the number of B ± mesons. The fit likelihood uses a Gaussian function 
to model the B^ signal. As in the previous measurements [9, 10], we use a first order 
polynomial to estimate the underlying combinatorial background. As shown by Fig. 1, the 
mass region above the B ± meson signal is not affected by partially reconstructed B ± decays 
and is quite well described by a straight line. We fit the data in the invariant mass range 
5.18 — 5.39 GeV/c 2 . The lower limit is chosen to avoid the region populated by partially 
reconstructed B ± decays. The width of the fitted mass range determines the statistical 
error of the background estimate. Since we have a much larger data set than previous 
measurements [9, 10], we can afford to fit the data in a smaller mass range in order to 
reduce the systematic uncertainty due to the background modeling. 

The average of the B ± mass values returned by the fits in the different px bins is 
5.2790 GeV/c 2 with a 0.5 MeV/c 2 RMS deviation, in agreement with the PDG value [25]. 
In the fit used to determine the number of B ± mesons, we fix the B ± mass value to 
5.279 GeV/c 2 [25]. The width of the Gaussian is a free fit parameter; the value of a 
returned by the fit increases from 12.0 ± 0.4 to 20.0 ± 0.4 MeV/c 2 from the first to last p T 
bin, in agreement with the simulation prediction. The fits are shown in Figs. 2 to 6. They 
return a signal of 2792 ± 186, 2373 ± 110, 1365 ± 66, 1390 ± 63, and 277 ± 44 B ± mesons in 
the five pt bins. 

We have investigated possible systematic uncertainties in the fit results. We have studied 
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FIG. 1: Invariant mass distribution of all candidates. The line represents a fit to the data 
using a first order polynomial plus a Gaussian function in order to estimate the background and 
the signal, respectively. 



the contribution of the B — ► J 'ftp n decay mode, the branching fraction of which is (4.9 ± 
0.6)% of that of the B — > J/ip K decay mode [25]. As shown in Ref. [26], the invariant mass 
distribution of Cabibbo-suppressed B decays, reconstructed assuming that pions are kaons, 
is shifted into the mass region 5.28 — 5.44 GeV/c 2 , which partially overlaps with that of 
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TABLE I: Detector acceptance, A, as a function of the px- The acceptance -4 cor r includes 
corrections evaluated using the data. The average < px > is the value at which the theoretical 
differential cross section [1] equals the integrated cross section in each momentum bin divided by 
the bin width. 

p T range (GeV/c) < p T > (GeV/c) A (%) A co „ (%) 

6 - 9 7.37 1.545 1.780 ± 0.045 

9 - 12 10.38 3.824 4.405 ±0.111 

12 - 15 13.39 5.966 6.872 ±0.173 

15 - 25 19.10 8.819 10.16 ±0.25 

> 25 12.516 14.42 ±0.36 



the B — > J /if) K decay mode. However, part of the Cabibbo-suppressed contribution is also 
used by the fit to predict the background under the B — > J /if) K signal with the effect of 
reducing its size. When adding the expected contribution of the Cabibbo-suppressed decays, 
the B — > J /if} K signal returned by the fit decreases by (1 ± 1)%. We have investigated 
other possible causes of systematic uncertainties in the B signal estimate. We have compared 
the results of our fit with those returned using an unbinned likelihood method. We have 
decreased the fitted mass range to 5.24 — 5.33 GeV/c 2 , and we have fitted the larger mass 
interval 5.18 — 5.60 GeV/c 2 . We have fitted the signal with two Gaussian functions in order 
to study detector resolution effects. The B signal returned by these fits does not vary by 
more than ±1.5%. Therefore, we attribute an overall ±2% systematic uncertainty to the fit 
results. 

A. Acceptances and Efficiencies 

The detector acceptance is calculated with a Monte Carlo simulation based upon a NLO 
calculation 2 . The B^ decay is modeled with the evtgen Monte Carlo program [27] that 
accounts for the J /if) longitudinal polarization [28]. The detector response to particles 
produced by B^ decays is modeled with the CDF II detector simulation that in turn is 
based on the geant Monte Carlo program [29]. The simulation includes the generation 
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FIG. 2: Invariant mass distribution of candidates with 6 < px < 9 GeV/c. The line represents 
the best fit to the data described in the text. 



of LI CMU trigger primitives. Simulated events are processed and selected with the same 
analysis code used for the data. The acceptances estimated using the simulation are listed 
in Table I. We use the data to verify the detector acceptance and efficiencies evaluated 
using the CDF II detector simulation. We adjust the simulation to match measurements 
in the data of: (1) the off-line COT track reconstruction efficiency; (2) the CMU detector 
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FIG. 3: Invariant mass distribution of B^ 1 candidates with 9 < px < 12 GeV/c. The line represents 
the best fit to the data described in the text. 



acceptance and efficiency; (3) the efficiency for finding LI CMU primitives; and (4) the 
efficiency of the LI, L2, and L3 triggers. 

In the simulation, the off-line COT track reconstruction efficiency is given by the fraction 
of tracks, which at generator level satisfy the pt and r\ selection cuts, that survive after 
selecting fully simulated events as the data. The COT track reconstruction efficiency is found 
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FIG. 4: Invariant mass distribution of candidates with 12 < px < 15 GeV/c. The line 
represents the best fit to the data described in the text. 



to be 0.998 ± 0.002. The same efficiency in the data is measured by embedding COT hits 
generated from simulated tracks into J /if} data. In Ref. [30], the COT track reconstruction 
efficiency in the data is measured to be 0.996 with a ~ 0.006 systematic accuracy 4 . We 

4 The efficiency measurement was performed in a subset of the data used for this analysis. Studies of 
independent data samples collected in the data taking period used for this analysis show that changes of 
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FIG. 5: Invariant mass distribution of candidates with 15 < px < 25 GeV/c. The line 
represents the best fit to the data described in the text. 



conclude that the efficiencies for reconstructing the /i + fi'K^ system in the data and the 
simulation are equal within a 2% systematic error. Kaon decay and interactions are modeled 
with the CDF II detector simulation. Because of the uncertainties of the detector materials 

the track reconstruction efficiency are appreciably smaller than the quoted systematic uncertainty [31]. 
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FIG. 6: Invariant mass distribution of candidates with px > 25 GeV/c. The line represents 
the best fit to the data described in the text. 



and the nuclear interaction cross sections, the kaon tracking efficiency has an additional 
0.3% uncertainty [32]. 

In the simulation, the fraction of CMU stubs generated by muon tracks with px > 
2 GeV/c and |r/| < 0.8 is 0.6439 ± 0.0004. In the data, this efficiency is measured by using 
J/ip — > fi + fi~ decays acquired with the /i-SVT trigger. We evaluate the invariant mass of 
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all pairs of a CMUP track and a track with displaced impact parameter, p T > 2 GeV/c, 
and \r)\ < 0.8. We fit the invariant mass distribution with a first order polynomial plus two 
Gaussian functions to extract the J/ip signal. From the number of J ftp mesons reconstructed 
using displaced tracks with or without a CMU stub (Fig. 7(a) and (b), respectively), we 
derive an efficiency of 0.6251 ± 0.0047. The integrated efficiency is evaluated after having 
weighted the pt and rj distributions of displaced tracks in the data to be equal to those of 
muons from B decays in the simulation. 

In the simulation, the efficiency for finding a CMU primitive (CMU stub matched by 
a XFT track) is 0.8369 ± 0.0004. This efficiency is measured in the data by using events 
acquired with the CMUPp^4 trigger. We combine the CMUP muon with all other CMU 
muons found in the event with and without a LI CMU primitive. We extract the number 
of J/ij} — > mesons by fitting the invariant mass distributions of all candidates with 

a first order polynomial plus two Gaussian functions. By comparing the fitted numbers of 
J ftp candidates with and without LI CMU primitive (Fig. 8(a) and (b), respectively) we 
derive an efficiency of 0.9276 ± 0.0005. The integrated efficiency is evaluated after having 
weighted the pt and rj distributions of the additional CMU muons to be equal to that of 
muons from B decays in the simulation. 

In the simulation, the efficiencies of the LI and L2 triggers are 0.9868 and 0.9939, respec- 
tively. By studying J/ip candidates acquired with the CMUPp r 4 trigger, the LI efficiency is 
measured to be 0.9879 ± 0.0009, and that of the L2 trigger 0.9948 ± 0.0001. The L3 trigger 
is not simulated. The L3 trigger efficiency is dominated by differences between the online 
and off-line reconstruction code efficiency 5 . The relative L3 efficiency for reconstructing a 
single muon identified by the off-line code has been measured to be 0.997 ± 0.002 [30]. The 
reconstruction efficiencies in the data and in the simulation are summarized in Table II. 

B. Results 

The differential cross section da/dp T is calculated as 



da(B+) 
dp T 



Ap T x C x ^l corr x BR 



N/2 



(1) 



5 Online algorithms are faster but less accurate than the off-line reconstruction code. 



20 





FIG. 8: Invariant mass distribution of a CMUP muon paired with all CMU muons in the event 
with (a) or without (b) a LI CMU primitive. 
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TABLE II: Summary of efficiencies for reconstructing B ± candidates in the data and the simulation. 
The last column indicates the corrections applied to the simulated acceptance and used to derive 
.A corr in Table I. 



Source 


Data 


Simulation 


Corr. 


COT tracking 


(0.996 ± 0.006) 3 


(0.998 ± 0.002) 3 


1.00 ± 0.02 


Kaon interaction 






1.000 ± 0.003 


CMU acc. and eff. 


(0.6251 ± 0.0047) 2 


(0.6439 ± 0.0004) 2 


0.942 ± 0.014 


Lf CMU primitives 


(0.9276 ± 0.0005) 2 


(0.8369 ± 0.0004) 2 


1.228 ± 0.002 


Lf eff. 


0.9879 ± 0.0009 


0.9868 


1.0011 ±0.0009 


L2 eff. 


0.9948 ± 0.0001 


0.9939 


1.0009 ± 0.0001 


L3 eff. 


(0.997 ±0.002) 2 


1 


0.994 ± 0.004 


Total 


0.328 ± 0.008 


0.283 ±0.002 


1.152 ±0.029 



where N is the number of B mesons determined from the likelihood fit to the invariant 
mass distribution of the J/ip K ± candidates in each bin. The factor 1/2 accounts for the 
fact that both B + and B~ mesons are used and assumes C invariance at production. The 
bin width Ap T and A corr , the geometric and kinematic acceptance that includes trigger and 
tracking efficiencies measured with the data, are listed in Table I. The integrated luminosity 
of the data set is C = 739 ± 44 pb _1 . The branching ratio BR = (5.98 ± 0.22) x 1CT 5 is 
derived from the branching fractions BR(B ± — > J/ip = (1.008 ± 0.035) x 10~ 3 and 
BR(J/4> -> fj+ fj,-) = (5.93 ± 0.06) x 10" 2 [25]. 

The measured B + differential cross section as a function of its transverse momentum is 
listed in Table III . The integrated cross section is 

a B+ fa > 6.0 GeV/c, \y\ < 1) = 2.78 ± 0.24 fib, (2) 

where the 8.8% error is the sum in quadrature of the 6% error on the integrated luminosity, 
the 3.6% uncertainty of the B + — > J/tp K + and J/ip — > fi + fi~ branching fractions, the 2.5% 
uncertainty of the acceptance calculation, the 2% systematic uncertainty of the fit, and the 
4.4% statistical error. 

For completeness, Figure 9 compares transverse momentum distributions in the data and 
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TABLE III: Observed differential cross section, da/dp? (nb/GeV/c), for B + mesons with rapidity 
\y\ < 1. Errors are the sum in quadrature of statistical errors (shown in parentheses) and systematic 
errors due to luminosity (6%), branching ratios (3.6%), acceptance (2.5%), and fitting procedure 
(2.0%). Systematic errors are not px dependent. The integrated cross section for px > 25 GeV/c 
is 21.7 ±3.7 nb. 



<pr> (GeV/c) 


Events 


Acceptance (%) 


da/dpT 


7.38 


2792 ± 186 


1.780 ± 0.045 


591.7 ± 59.0 (39.3 stat. ) 


10.38 


2373 ± 110 


4.405 ±0.111 


203.2 ± 17.8 (9.4 stat.) 


13.39 


1365 ± 66 


6.872 ±0.173 


74.9 ± 6.6 (3.6 stat.) 


19.10 


1390 ± 63 


10.16 ±0.25 


15.5 ± 1.3 (0.7 stat.) 


> 25 


277 ± 44 


14.42 ± 0.36 





in the simulation, based on the NLO QCD prediction 2 , that has been used to evaluate the 
detector acceptance. Data and simulation are normalized to the same number of events. 
Each distribution is constructed using J/i[)K ± candidates with invariant mass in the range 
5.255 — 5.315 GeV/c 2 (region #1). The background contribution is subtracted using candi- 
dates in the mass range 5.18 — 5.24 and 5.33 — 5.425 GeV/c 2 . The background normalization 
is the number of events in region #1 minus the number of B ± candidates determined by 
the fit listed in Table III. One notes the fair agreement between data and untuned QCD 
prediction 6 . 

V. CONCLUSIONS 

We use the exclusive decay B ± — > J/i/j K ± to measure the B + production cross section 
in pp collisions at y/s = 1960 GeV. The measurement is based on a sample of 8197 ±239 B ± 
mesons selected from 739 pb _1 of data collected with the CDF II detector at the Fermilab 

6 The pt distributions of the B ± and J/ip mesons in the data are slightly softer than those of the simulation; 
this difference is not relevant for the result of the study because the kinematic acceptance has been 
evaluated for each p^-hm and the calibration of the simulated acceptance using the data do not depend 
on the muon and kaon transverse momenta. 
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FIG. 9: Transverse momentum distributions in the data (•) and simulation (solid histogram). The 
simulation is normalized to the — » /i + n~ signal observed in the data (see text). 



Tevatron collider. The B + production cross section is measured to be 

<j b +{pt > 6.0 GeV/c, \y\ < 1) = (2.78 ± 0.24) ^b. 

To compare with other Tevatron measurements, we choose as a theoretical benchmark 
the NLO QCD prediction [1] that uses a 6-quark mass of m& = 4.75 GeV/c 2 , renormalization 



and factorization scales [Ir = jXp = Jpt + m h the MRSD [6] fit to the parton distribution 
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functions (PDF), a fragmentation fraction /„ = 0.375, and a fragmentation model based on 
the Peterson fragmentation function with the e parameter set to 0.006. The ratio of the 
present measurement to this theoretical prediction is 2.80 ± 0.24. Previous measurements 
of the single 6-quark cross section based on the detection of J/ip mesons yield the following 
ratios to the same theoretical prediction: 2.9 ± 0.67 [9], 4.0 ± 0.6 [10], 4.0 ± 0.4 [33], and 
3.14 ± 0.28 [30]. In contrast, all CDF and DJZ) measurements of the single b production 
cross section that are based upon detection of a lepton from 6-quark decays [11-15] yield a 
smaller average ratio to the same theoretical prediction (2.2 with a 0.2 RMS deviation [5]). 
As shown in Fig 10, our measurement agrees with the value inferred from the J/ip inclusive 
cross section [30] [&b+{Pt > 6.0 GeV/c, \y\ < 1) = 2.4 ± 0.4 fib] and is within the range 
of values predicted by the FONLL QCD calculation [2, 34] that uses f u = 0.389 [25] and 
the CTEQ6M fits to the parton distribution functions [35] (2.1 fib with a ~ 30% theoretical 
uncertainty [4]). 
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FIG. 10: Measurements of the B + differential cross section (\y B | < 1) at the Tevatron are 
compared to the NLO and FONLL theoretical predictions (see text). The result of this experiment 
(•) is shown together with those of (A) Ref. [30] and (o) Ref. [9]; the result of Ref. [9] has been 
increased by 10% to account for the expected increase of the cross section from y/s = 1.8 to 1.96 
TeV [1]. 
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